Zeros of some bi-orthogonal polynomials. 
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Abstract. Ercolani and McLaughlin have recently shown that the zeros of the bi- 
orthogonal polynomials with the weight w{x, y) = exp[—{Vi{x) +V2{y) + 2cxy)/2], relevant 
to a model of two coupled hermitian matrices, are real and simple. We show that their 
argument applies to the more general case of the weight {wi*W2*-..*Wj){x, y), a convolution 
of several weights of the same form. This general case is relevant to a model of several 
hermitian matrices coupled in a chain. Their argument also works for the weight W{x, y) ~ 
Q-^-y I (^x + y), < x,y < CO, and for a convolution of several such weights. 
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1. Introduction. For a weight function w{x,y) such that all the moments 



y)x^y-' dxdy 



exist and 

Dr, := det[M,,,],,,=o,i,...,n ^ (L2) 

for all n > 0, unique monic polynomials Pnix) and Qnix) of degree n exist satisfying the 
bi-orthogonality relations (a polynomial is called monic when the coefficient of the highest 
degree is one) 

w{x,y)pn{x)qm{y)dxdy = hnSmn- (1-3) 



Just like the orthogonal polynomials they can be expressed as determinants, e.g. 



D 



det 



n-l 



rMo,o 

Mi,o 



Mo,n-l 1 

Mi,n-i X 



and have integral representations, e.g. 

» n 

« y A„(x)A,.(y) n(- wMV!/. 



An(x):= n (^^-^0, A„(y):= J] {y, 



Vi] 



;i-4) 



(1-5) 
(1-6) 



l<i<j <n 



l<i<j<n 
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From limited numerical evidence for the weights 

(i) w{x, y) = sm{nxy), < x, y < 1; 

(ii) w{x, y) = \x-y\, -l<x,y<l; 

(iii) w{x, y) = [l/{x + y)] exp[-x - y], < x, y < oo; 

(iv) w{xj y) = exp(— a;^ ~ ~ cxy)^ — oo < x,y < oo, < c < 2; 

one might think that the zeros of the bi-orthogonal polynomials are real, simple, lie re- 
spectively in the x or y-support of w{x,y), interlace for successive n, ... 

Alas, this is not true in general as seen by the following example due to P. Deligne. 
If one takes 



w{x,y) = u{x,y) + v{x,y), 

S{x-y), -l<x,y<l, 
0, otherwise, 

vix,y) = ^[5{x - l)5{y + 2) + Six + - 2)]. 



u{x,y) 



(1.7) 
(1.8) 

(1.9) 



Then the zeros of P3{x) and q3{x) are complex. 

However, N.M. Ercolani and K.T.-R. Mclaughlin have recently [1] shown that with 
the weight function 



wi{x,y) = exp 



■\vi{x) - ^V2{y) - cixy 



(1.10) 



(— oo < x,y < oo), Vi and V2 polynomials of positive even degree, c a small non-zero 
real constant, all the zeros of the bi-orthogonal polynomials Pnix) and qnix) are real and 
simple. 

In this brief note we will show that their argument works for the following general case 
encountered for random hermitian matrices coupled in a linear chain. Let Vj{x), ^ ^ j ^ P, 
be polynomials of positive even degree and cj, 1 < j < p, be small real constants, none 
of them being zero ("small" so that all the moments Mij defined below, eq.(1.13), exist.) 
Further let 

-Ti^kix) - - Vk+i{y) - Ckxy 



Wk{x,y) := exp 



(1.11) 



2-^/2 
Moreover, assume that for alH, j > 



^iJ '■= J x\wi *W2* ... * Wp-i){x, y)y^dxdy 



(1.13) 



exist. 



Theorem. Then monic polynomials pj (x) and qj (x) can be uniquely defined by 

J pj{x){wi *W2* ... * Wp-i){x,y)qk{y)dxdy = hjSjk (1-14) 
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and all the zeros of pj (x) and of qj (x) are real and simple. 

The same argument works for any weight W{x, y) such that det[VF(xi, |/j)]i,j=i,...,n > 
for xi < a:2 < ... < x^, Hi < y2 < ••• < Vn and moments M^j = f W{x,y)x^y^dxdy exist 
for all z,j > 0. For example, if W{x,y) = [l/{x + y)] exp[— a; — y], < x,y < oo, then 
monic polynomials pj {x) can be uniquely defined by 



/ pj{x)W{x,y)pk{y)dxdy ^ hjdjk (1.15) 

^0 



(here W{x,y) is symmetric in x and y so that Pj{x) = qj{x)) and all the zeros of Pj{x) are 
real, simple and non-negative. 

2. Results and proofs. Here we essentially follow section 3 of reference [1]. With 
any monic polynomials Pj{x) and qj{x) of degree j, let us write 

P,j{x) := pj{x) (2.1) 

Piji^) '■= J Pj(C)(«^l *W2* ...*Wi_i){^,x)d^ 

■■= JpjmLi{^,^)d^, l<i<p (2.2) 
Qp,j{x) •■= qj{x) (2.3) 
Qiji^) — J {wi*Wi+i * ...*Wp_i){x,^)qj{^)d^ 

:= J UR,{x,C)qj{m l<i<P (2.4) 

Lemma 1. For xi < X2 < ■■■ < x^, yi < y2 < ••• < Vn, 

det K(a;j,yfc)]^- > 0. (2.5) 

This is essentially eq. (40) of reference [1]. This can also be seen as follows. Let X — [xidij] 
and Y = [yiSij] be two nxn diagonal matrices with diagonal elements xi, Xn and j/i, 
yn respectively. Then the integral of exp[— ctr t/Xt/~-^y] over the nxn unitary matrices 
U is given by [2] 

det [exp(— cxi-y,)! • ■ , 

K — ^ , ^y^y=i,...n ^2.6) 

A,(x)A,(y) ^ ' 

where is a positive constant depending on c and n. Hence 

\Y,{Vi{x^) + Vi+^{y,)) 



exp 



2 



,^,-e.t..x.-^.,^detK(x,,.)] 

A„(x)A„(y) ^ ' 



The left hand side is evidently positive while on the right hand side the denominator is 
positive since xi < X2 < ■■■ < Xn and yi < y2 < ■■■ < yn- From this eq. (2.5) follows. 
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Lemma 2. For xi < X2 < ... < Xn, Hi < 1/2 < ■■■ < Vn, 

det [{wi, * Wi^ * ... * Wi^){xj, yk)]j^k^i^,„^n > 
Proof. Binet-Cauchy formula tells us that [3] 

det [{wi, * Wi^){xj,yk)]^^k=i,...,n 

is equal to 



(2.8) 



^i<6<...<^„ 



det [wi, {xj, 6)]j,fe=i,...,n • det [wi^i^j, yfe)]j-,fe=i,...,„ d^i-d^n (2.9) 



By lemma 1 the integrand is every where positive, so lemma 2 is proved for the case £ = 2. 
The proof is now completed by induction on £, using again the Binet-Cauchy formula. 

Lemma 3. For any monic polynomial Pj{x) of degree j, Pij{x), 1 < i < p, has at 
most j distinct real zeros. Similarly, for any monic polynomial qj{x) of degree j, Qij{x), 
1 < i < p, has at most j distinct real zeros. 

Proof. Let, if possible, zi < Z2 < ... < Zm, rn > j,he the distinct real zeros of Pij{x). 
Since 

j 

Pi,j{x) = J2^kTi,k{x), (2.10) 

k=0 



with 



Pidi^i)^^^ £ = l,2,...,m, m>j 



(2.11) 
(2.12) 



implies that 



= det 



Tifi{zj+i) Ti^i{zj+i) ... Tij{zj+i) 



= / det 



ULi{^l,Zi) ^2ULi{^2,Zi) ... ^j+iULi{^j+i, Zi) 
_ULi{^l, Zj + i) C2ULi{C2,Zj+i) ... Cj^iULi{ij + i,Zj+i) _ 



d^i...d^j+i 

(2.13) 



or 



J det[UUCk,Zi)]„^,^,_j+,.det [^t']J,^,^,_J+,dC,...d^,+^ = (2.14) 



in contradiction to lemma 2. Thus m can not be greater than j. 
The proof for Qij{x) is similar. 



Lemma 4. Let the real constants ci, Cp-i, none of them being zero, be such that 



exist for all i, j > 0. Then 



Dn := det[M,,,],,,=i,...,^ ^ 



(2.15) 
(2.16) 



for any n > 0. 



and 



or 



Proof. Let, if possible, — for some n. Then Yl^=o ^ij^lj — ^-i zero. 



/n 
x'-ULp{x,y)^qjy^dxdy = 0, i = 0, 1, 

/n 
Piix)ULp{x, y) ^ qjy^dxdy = 



(2.17) 



(2.18) 



for any polynomial Pi(a;) of degree i <n. But 



/n 
ULp{x,y)^qjy^dy 

has at most n distinct real zeros (lemma 3). So one can choose Pi{x) such that 

/n 
ULp{x,y)^qjy^dy > 
j=0 



(2.19) 



(2.20) 



in contradiction to eq. (2.18). So Dn ^ and bi-orthogonal polynomials (x) , qj{x) exist, 
see eqs. (1.4), (1.5). 

Lemma 5. Let Pj{x), qj{x) be the bi-orthogonal polynomials, eq. (1-14); or with the 
definitions (2.1)-(2.4) 



Pi,j{x)Qi,k{x)dx ^ hj5jk, '^<i<P 



(2.21) 



Then Pij{x) has at least j real distinct zeros of odd multiplicity. So does have Qij{x). 

Proof. Let, if possible, zi < Z2 < ... < Zm, rn < j, be the only real zeros of Pij{x) of 
odd multiplicity. Set 



R{x) = det 



Qi,oix) Qt,i{^) 
Qi,o{zi) Qi,i{zi) 



Qi,oiZm) Qi,l{Zm) 

m 

Um{x,0j2''k^'dC 



Qi,m ('^) 



(2.22) 



(2.23) 



with some constants depending on zi, Zm- 
Since m < j, the bi-orthogonahty gives 



Pij{x)R{x)dx = 0. 



(2.24) 



However, R{x) can also be written as 



= / det 



det 



■ URi{x,^o) 


URi{x, 


URi{zi,^o) 


URi{zi,^l)^l 


-URi{Zm,^o) 


URi{Zm, ^1)^1 


' URi{x,^o) 


URi{x,^l) 


URi{zi,^o) 


URi{zi,^l) 


-URi{Zm,^o) 


URi{Zm,^l) 



URi{x,$,m)C 
URi{zi,^rn)i 



m -, 

m 
m 
m 



Uri {Zm ) Cm)Cm - 
URi{x,^rn) 
URi{zi,^rn) 



(m + 1)! 



det 



Uri {Zqji, Cm) - 

URi{x,^o) URi{x,Ci) ... URi{x,^r 
URi{zi,^o) URi{zi,^i) . 



det 



-URi{Zm,^o) URi{Zm-,il) 

URi{x,io) UR^{x,ii) 

URi{zi,Co) URi{zi,^l) 



eo<^i<---<^r, 



URi{zi,^rn) 

m J 

URi{x,im) 
URi{zi,^rn) 



n iCs - Cr)dCod^i...d^, 



0<r<s<m 



URi{Zm,^o) URi{Zm,Cl) ■■■ URi{Zm,Cm) 



n i^s-^r)d^0d^l...d^r, 

(2.25) 



0<r<s<m 



Thus R{x) is represented as an integral whose integrand has a fixed sign determined by 
the relative ordering of the numbers x, Zi, Z2, Zm (lemma 2). It thus follows that R{x) 
changes sign when x passes through any of the points Zk, k = 1, m and at no other value 
of X. In other words, zi, Zm are the only real zeros of R{x) having an odd multiplicity. 
And therefore Pij{x)R{x) has a constant sign, so that 



Pij{x)R{x)dx ^ 



(2.26) 



in contradiction to (2.24). 

The proof for Qij{x) is similar. 

As a consequence we have the integral representations of Pi,j{x) for z > 1 and of 
Qij (x) for i < p involving their respective zeros 



Pij{x) cx / det 



ULi{^0,x) UlMi,x) ... ULi{^j,x) 

ULi{^0,Zi) ULi{^i,Zi) ... ULi{(,j,Zi) 

ULi{^0,Zj) ULi{^i,Zj) ... ULi{^j,Zj). 
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n {Cs-Cr)dCodCi...dCj 

0<r<s<j 



Qij{x) oc / det 



URi{zi,^o) URi{zi,^i) ... URi{zi,^j) 

-URi{Zj,^0) URi{Zj,^l) ... URi{Zj,^j) 



(2.27) 



0<r<s<j 

(2.28) 



Lemmas 3 and 5 tell us that if Pj{x) and qj{x) are bi-orthogonal polynomials satis- 
fying eq. (1.14), then Pi_j{x) and Qij{x) each have exactly j distinct real zeros of odd 
multiplicity. In particular, the zeros of the bi-orthogonal polynomials Pj{x) = Pij{x) and 
qj{x) = Qpj{x) are real and simple. 

With a little more effort one can perhaps show that all the real zeros of Pij{x) and 
of Qij{x) are simple. Other zeros, if any, must be complex. Whether the zeros of Pj{x) 
{qj{x)) interlace for successive j, remains an open question. 

3. Bi-orthogonal polynomials v^ith another weight. 

For the weight W{x, y) — [l/{x + y)] exp[—x — y], < x,y < oo, one can say as follows. 
Lemma 1'. One has [4] 



det[W{xj,yk)]j,k=i,...,n = exp 



A„(x)A,(y) II (xj+yk)-' (3.1) 
j,k=i 



which is evidently positive for < xi < X2 < ... < x^, < yi < y2 < ■■■ < y-n- 

Lemma 3'. For any monic polynomial pj(x) of degree j, Pj(a;) := W{x,y)pj{y)dy 
has at most j distinct real non-negative zeros. 

In the proof of lemma 3, replace eqs. (2.10)-(2.14) by 



= det 



-f 

Jo 



Toizi) 



k=0 

oo 



Tkix)= / eW{^,x)dC 
Jo 

Pj{zi)=0, £=1,2, ...,m, m>j 
Ti(zi) ... T,(zi) 



To{zj+i) Ti{zj+i) ... Tj{zj+i) 
W{Ci,zr) ^2W{^2,zi) 



det 



ej+iW^(e.-+i,^i) 



W{Cl,Z, + i) C2W{C2,Zj + l) ... e+iW{Cj + i,Zj + i) 



(3.2) 

(3.3) 
(3.4) 



d^i-.-d^j+i 



J det [W{^k, Z£)]k,£=i,...,j+i d^i-d^j+i 



(3.5) 



ur 

/ det[W-(a,^^)]fc,,=i,...,,+i.det [^i-']k,i=,_j+,d^i-d^j+i = 

in contradiction to lemma 2'. Thus m can not be greater than j. 
Lemma 4'. With 

/■oo 

Mjj := / x''W{x,y)y^dxdy, 
Jo 

Dn := det[Mij]jj=o,i,...,n 7^ 

for any n > 0. 

In the proof of lemma 4 replace everywhere / Ulp{x, y)... by W{x^y). 
Lemma 5'. Let Pj{x) be the (bi-orthogonal) polynomials satisfying 



(3.6) 

(3.7) 
(3.8) 



f 

Jo 



W{x,y)pj{x)pk{y)dxdy = hjSjk- 



(3.9) 



Then Pj{x) := W{x, y)pj{y)dy andpj(a;) each have at least j distinct real non-negative 
zeros of odd multiplicity. 

Let, if possible, Q < zi < Z2 < ... < Zm tu < j ,he the only real non- negative zeros of 
Pj{x) of odd multiplicity. Set R{x) = YYjLii^ ~ ^j)- Then as m < j, one has 



f 

Jo 



Pj{x)R{x)dx = 



(3.10) 



But Pj{x) and R{x) change sign simultaneously as x passes through the values zi, 
Zm and at no other real positive value. So the product Pj{x)R{x) never changes sign, in 
contradiction to (3.10). Therefore Pj{x) has at least j distinct real non-negative zeros of 
odd multiplicity. 

To prove that Pj{x) has at least j distinct real non-negative zeros let if possible, 
< zi < Z2 < ... < Zm, Tn < j, be the only such zeros. Set 



R{x) = det 



Jo 



Po{x) Pi{x) ... PM 

Po{zi) Pl{zi) ... Pm{zi) 

iPo 

m 

w{x, e) J2 ^k^d^ 

k=0 



(3.11) 



with some constants depending on zi, Zm. 
Since m < j, the bi-orthogonality gives 



Pj{x)R{x)dx = 0. 



(3.12) 



But 



R{x) oc / det 



W{x,^o) W{x,Ci) ... W{x,U) 



n {^s-Cr)d^od^i...dC. 



0<r<a<m 



(3.13) 



which says that Zi, Zm are the only distinct real no n- negative zeros of R{x) and therefore 
Pj{x)R{x) has a constant sign, in contradiction to (3.12). 

Conclusion. We have shown with the arguments of Ercolani and McLaughlin that 
if the weight w{x,y) is such that dct[w{xi,yj)]ij^i^^^^^n > for xi < X2 < ... < Xn, 
Hi < y2 < ■■■ < Un and moments J w{x,y)x'^y^ dxdy exist for all i,j > 0, then bi-orthogonal 
polynomials exist and their zeros are real, simple and lie in the respective supports of 
the weight w{x,y). The same is true for a weight which is a convolution of several such 
weights. 
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